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ABSTRACT 



In this thesis, large deflection plate compatibility 
and equilibrium differential equations are derived using finite 
deformation theory including the non-linear five constant stress- 
strain relationship. In addition, the assumptions of classical 
large deflection plate problems under classical elasticity 
theory are discussed in connection with the more general 
approach proposed by this thesis. 
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INTRODUCTION 



The approach to the large deflection plate problem pro- 
posed by this thesis, as suggested by the title, is one utilizing 
non-linear elasticity theory, or theory dealing with finite defor- 
mation of elastic solids. In reference [l] Francis D. Murnaghan 
provides an excellent discussion of the theory, and derives the 
basic stress-strain relationship used in this thesis (see appendix 
B) . The use of this basic non-linear stress-strain relationship, 
sometimes called Murnaghan' s law |_2], places the scope of this 
thesis more appropriately into the field of rheology rather than 
structures per se (see Figure l). 

In actuality, predicting how materials will behave in 
response to forces lies in the domain of rheology - literally 
the study of flow. The Society of Rheology uses the words of 
Heraclitus, ITO&vra (everything flows), as its motto, and 

books on rheology quote this philosophical conviction in their 
introductory pages. [3]. Because rheology deals with the flow 
of matter under the action of forces, by an obvious extension 
it also includes all deformations of materials by forces. Thus, 

Li] References are listed beginning on page 66. 

L 2 ] Novozhilov states "The elastic law corresponding to it (the 
five constant stress-strain theory) is ordinarily called 
Murnaghan 's law, although it was actually first proposed 
much earlier by Voigt in 1893. The first attempt to examine 
the stress-strain relation in a form different from Hooke's 
law was made by Bulffinger in a paper published in the works 
of the Russian Academy of Sciences in 1729" pp. 127. 
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as shown in Figure 1 , rheological theories have as their model 
the classical theories of elasticity and hydrodynamics in which 
it is assumed that a material, whether solid or liquid, deforms 
linearly in response to a force. From the classical models have 
grown more and more general theories, first by abandoning classical 
distinctions between a solid and a liquid, and then by abandoning 
the requirement that the response be linear. The most general 
theory, as shown in Figure (l), is the theory of nonlinear visco- 
elasticity which includes all other theories as special cases. 

A theory of non linear visco-inelasticity exists (shown in Figure 
1 dashed) but is of historical interest only [3]. 




FIGURE 1 

(From "Deformation and Flow" by Charles J. Lynch as appeared in 
International Science and Technology January 1966) 
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Since the end of World War II there has been a consid- 



erable change in the approach to rheology and to the formulation 
of rheological theories. The origin of this change is two-fold. 
First, certain flow phenomena, observed first in flame-thrower 
fuels and later in a wide variety of other fuels, were in direct 
opposition to classical hydrodynamic theory. The other factor, 
which was perhaps of greater significance in its contribution 
to the changing approach to rheology, at least to those primarily 
interested in structures, was the development of a continuum 
theory for elastic materials. 

Although the theory of finite deformation of an elastic 
solid applies equally to all elastic materials, most of the 
research accomplished to date has been in the rubber industry. 

The theory as applied to metallic structures has only recently 
become of practical interest with the increased emphasis on 
limit design, and the more stringent requirements of high per- 
formance aircraft and missiles, space craft, and deep submergence 
pressure hulls. In the majority of experimental research con- 
ducted so far, rubber is generally used because it has a large 
elastic range and its nonlinearities are therefore more easily 
measured . 

Vulcanized rubber differs from other elastic materials in 
the extent of its elastic deformation. A rubber rod or strip 
may be stretched to four or five times its initial length without 
suffering permanent deformation. When rubber is subjected to 
shear deformation and the corresponding shearing force is measured, 
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the relation between the shear strain and the shearing force 
is found to be approximately linear, just as with most elastic 
solids. On the other hand, when rubber is stretched by tensile 
forces, the relation between the tensile force and the elonga- 
tion shows considerable departure from linear behavior. 

While this may at first sight appear surprising, it can 
in fact be demonstrated mathematically that a material having 
a linear shear relationship cannot possibly have a linear tensile 
relationship. From this we can conclude that a true Hookean 
solid is only a convenient fiction - an elastic material cannot 
possibly have a linear relationship in both tension and shear as 
assumed in classical elasticity theory. The classical theory is 
only valid because the departure from linearity does not show up 
when the deformations are small - even a curved line appears 
straight if one is only concerned with a little piece of it. [3j 

In finite elasticity theory, the material properties are 
characterized by an expression for the energy of deformation per 
unit initial volume, ( M) ), [l] (see appendix B). This strain 
energy function and the associated material constants have been 
determined by experiment for various vulcanized rubbers and the 
actual manner in which stored energy depends on the deformation 
has been determined from them. Using this experimentally deter- 
mined stored energy, the results of other experiments with vulcan 
ized rubber have been predicted with considerable accuracy, pro- 
viding an excellent verification of the theory [3]. Although 
the material constants of interest to the structural engineer 
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(modulus of elasticity and shear modulus plus the three addi- 
tional third order elastic constants shown in appendix B) have 
been evaluated for several structural materials as shown in 
Table I, the numerical values show a wide variation between 
investigators with a wide spread for some materials and much 
more work is required. 



TABLE I 

THIRD ORDER ELASTIC CONSTANTS ( i numbers are the range of values) 



Reference and 
Material 


E psi X 10" 6 

Relative Value 


psi X 10 -6 


m 

psi X 10 -6 


n c 

psi X 10 0 


Hughes and Kelly 
Armco Iron 


27-5 


- 34.8 

± 6.5 


- 10.3 

± 7-0 


+ 110 
± 110 


Hughes and Kelly 
By rex 


2.0 


4 1.4 
± 4.0 


+ 9.2 
t 5.0 


4 42 
± 35 


■ ,j 11 1 \ 

Seeger and Buck 
Copper 


17.0 


- 22.6 
± 10 


- gg.o 

± 1.3 


- 225 
± 3 


Seeger and Buck 
Iron 


29-5 


- 24.3 
± 6 


- 110 
± 1.5 


- 215 
t 1.5 


Smith 

.6 Carbon Steel 


21.1 


- 46 


- 60 


- 67 


Smith 

Austenitic Steel 


18.3 


- 53.5 


- 75.2 


- >40.0 j 

I 

| 


Crecraf t 
Nickel-Steel 


17.4 


- 4.6 


- 59 


- 73-0 


Rollins etal. 

601 6-T6 Aluminum 








- 31.2 


Borg - 

Calculated values 
based on 
V = .25 


30.0 




-45 


- 34.2 



Data as reported in reference |_6] 
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Finite elasticity theory has successfully predicted a 
number of effects that are not predicted by classical theory. 

For example, experimental observations of a rod of circular 
cross section subjected to simple torsion made by Polynting 
in 1909 with steel wires and in 1913 with vulcanized rubber 
rod showed that the rod will not only twist as predicted by 
classical theory, but will also elongate. As is often the 
case when knowledge available in the literature remains unused 
for some time, these results seem to have been forgotten until 
recently when the application of finite deformation theory 
predicted this result [3]. 

That Murnaghan's finite deformation theory [l] should 

> 

be applied to plates and shells in the large deflection region 
has been proposed verbally by S. F. Borg for some time. In refer- 
ence [4], Borg, Hoppe and Kopchinski, and in references [5], [6], 
[ 7 ] and [8] Borg discusses large deflection of plates, including 
the derivation of plate compatibility equations from a finite 
deformation theory approach. This thesis is a continuation and 
extension of Borg's work as applied to plates. 
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GENERAL DISCUSSION OF NON LINEAR ELASTICITY 



Non linearity is introduced into the theory of elasticity 
dealing with isotropic materials in three ways. It can be intro- 
duced through the generalized strain tensor or strain displacement 
relation (appendix A), through the stress-strain relations (appen- 
dix B) or through the equations of equilibrium of a volume element 
of the body (appendix C). 

In the equations of equilibrium and the strain tensor 
relations, the retention of non-linear terms is conditioned by 
geometric considerations, or the necessity of considering the 
angles of rotation in determining dimensional changes of line 
elements and in formulating the conditions of equilibrium of a 
volume element. However, in the stress-strain relation, non- 
linear terms appear if the strain exceeds in magnitude, physical 
constants characteristic of the material and referred to as the 
limits of proportionality. It should be noted that the above 
discussion does not presuppose a linear elasticity origin. 
Actually, linear theory is developed by successively ignoring 
those terms or quantities which give non-linear characteristics 
to the resultant equations. For many materials, e.g. aluminum, 
the limit of proportionality can be quite low or non existant, 
and a finite deformation, continuum, theory is required where 
deformations are sufficiently large as to invalidate linear 
elasticity theory (Hooke's law). 

Thus, two types of non-linearity must be considered in 
large deflection problems, geometrical and physical. Since, in 
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general, smallness of angles of rotation of the body as a whole 
does not imply the smallness of elongation or shear strain of an 
infinitesimal element of volume of the body (and conversely), the 
geometrical and physical non-linearity can be regarded as inde- 
pendent of each other. This independence is useful for purposes 
of classifying problem types and for purposes of making reasonable 
assumptions to simplify the problem to enable an engineering 
solution. 

As a result, elasticity problems can be classed as one of 
four basic types [ 2 ]: 

(1 ) Elasticity Problems having both physical and geometrical 
linearity . 

In problems of this type, the effect of the angles of 
rotation are of the same order of magnitude as the elongations 
and shear strains, while the elongations do not exceed the limit 
of proportionality of the material. An example of this type of 
problem is the simple tension test when the stresses are main- 
tained below the proportional limit. 

(2) Elasticity Problems which are physically non linear 
but geometrically linear . 

For these problems, angles of rotation can be neglected 
in projecting the forces which act on a volume element and in 
determining strains. However, 'the elongations exceed the pro- 
portional limit and require a non-linear stress-strain relation. 
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The simple tension test becomes a problem of this type when the 
stresses in the rod exceed the proportional limit. 



( 3 ) Elasticity Problems linear physically , but non - 
linear geometrically . 

In problems of this type, the angles of rotation are 
essentially large but the strains do not exceed the limit of 
proportionality. An example of this type is the bending of a 
thin steel strip, or the buckling of a slender column within 
the elastic range where the original shape and position is re- 
gained after removal of the load, i.e. no permanent deformation. 

(4) Elasticity Problems non linear both physically and 
geometrically . 

In problems of the fourth type, the strains exceed the 
limit of proportionality and the angles of rotation are so large 
that it is necessary to retain non-linear terms both in the stress- 
strain relation, and in the equations of equilibrium. Plate and 
shell problems become examples of this type if the deformation 
is large and the stresses exceed the limit of proportionality. 
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LISCUSSION OF PLATE PROBLEM 



A plate may be defined as a three-dimensional deformable 
body whose thickness is small compared to its least lateral 
dimension. Because the thickness is small compared to the 
other linear dimensions, classical elasticity theory treats a 
relatively thin plate in an essentially two dimensional manner. 

The linear small-deflection theory of plates, developed 
by Lagrange (1811) is based on the following assumptions 1.9]: 

1 . points which lie on a normal to the mid-plane of the 

undeflected plate lie on a normal to the raid-plane of the deflected 
plate ( $| » jj « 0 ) s 

2. the stresses normal to the mid-plane of the plate, 

arising from the applied loading, are negligible in comparison 
with the stresses in the plane of the plate ( ) 6^ ) ; 

3. the slope of the deflected plate in any direction is 

small so that its square may be neglected in comparison with 
unity (curvatures ~ ) ; 

r» a* 1 5 ’ 

4. the mid-plane of the plate is a neutral plane: i.e. 

any mid-plane stresses or membrane stresses arising from the 
deflection of the plate into a non-developable surface may be 
ignored . 

In classical large deflection theory where deflection is 
of the order of magnitude of the plate thickness, and greater, 
the first three assumptions above for small deflection theory are 
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retained. However, the fourth assumption is no longer valid with 
large deflection, and middle surface stresses must be considered. 

As a result of the large deflection, the problem becomes 
non linear. The question is, does the problem become one of 
geometric non linearity as assumed in classical theory, or is it 
one of both physical and geometrical non-linearity as proposed 
by this thesis? Physical non linearity is currently introduced 
in the inelastic buckling problem as a variation in Poisson's 
ratio and variation of the Elastic Modulus in the so-called tan- 
gent/secant modulus theory. Gerard, [lO] provides a good "state 
of the art" summary of this theory. However, it is highly possible 
that the correct approach lies in the determination of the addi- 
tional elastic constants introduced by finite deformation theory. 

Assuming, for the moment, that the large deflection plate 
problem is physically linear, but geometrically non linear, there 
is some error introduced by the classical differential equations 
formulated by vonKarman [ 1 1 ] due to his treatment of the problem 
in a two dimensional manner. This error is negligible for very 
thin plates but becomes significant for plates of finite thickness. 
To illustrate this error, the classical large deflection assumptions 
will be applied to the generalized large deflection strain tensor 
(appendix A) . 

The assumption that points which lie on a normal to the 
mid-plane of the undeflected plate lie on a normal to the mid- 
plane of the deflected plate require that = 0 • 



Thus the large deflection strain tensor derived in appendix A 
reduces to: 




or 



j_ r^v , ^ , i/^^ 

z L di j j 



In the von Karman formulation; 



— \ Z /^i \ 2 " 

* + > > (fyj , <^X <5^ j 



&V /t^Vx 1 - , /£v\ 2 - 

> (axj and C^J ) 



are neglected, which seems reasonable, but even without the 
above terms, there are terms remaining in the expressions for 
^X2, ¥«J* and which are still significant. 

If we assume that the elements of the generalized strain 
tensor, as applied to a plate of finite thickness, are of the 
order of magnitude shown below, i.e. 



(l) Terms and are of the order €. and are 

<?X c5ij 

of the same order of magnitude as the derivatives of the 
shear forces Q x and Qy . 
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and 



(2) 



Terms 



a* 

a/ > 5^ > 9 >% } ^ 



/ <^ w \ 

%>z. 1 \5*/ J / 



2 . 

$tt ^ are of the order of magnitude 6 

d* 



(3) Terms || , >(lj) > dij) , 35 ' tx ^ > (J^ ) 

) Z and^^dH^ 2, are of the order of magnitude 



and can be neglected. 

Retaining terms up to the order of magnitude £ the generalized 

large deflection plate strain tensor can be reduced to 



( 2 ) 






5* + IteiJ'f 

Z\dA/ 


1 /^u . 

2^ 55 ^ 


1 £*/ 

2 a* 


-L/5u + & 

1 «ix 




i_ ^ w 
2. <? ^ 


1 

V 2 d* 


A, 

Z d ij 


< ^ 2: 



\ 



If the Poisson’s effect across the thickness of the plate is 
ignored, and ^ is taken to be zero, the above strain tensor 
becomes identical to that obtained by Borg, Hoppe, and Kopchinski L 4 J* 
In the differential equations formulated by von Karman L 1 I ] , 
the strain elements are assumed to be 



(3) 


= 






\ z 

SJ5x>> 


(4) 


e 3 = 


^3 


+■ T( 


K ^) 


(?) 


it 

3~> 

>2 


<$U 




J, 
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By comparing the von Karman elements with the elements of the 
three dimensional plate strain tensor, it is shown that von Karman 

i 

retains terms of the order £. and yet neglects terms of the 
order €. which are introduced by the shear strain elements Y* st- 
and Y^z , and also ignores the Poisson's effect across the 
thickness of the plate . Thus, even with the assumption 

that the problem is linear physically and non-linear geometrically, 
a feeling for the magnitude of the error resulting from the two- 
dimensional approach can be achieved. This points out a significant 
value of the continuum mechanics approach to the problem. 

At this point, it should be noted tliet the large deflection 
strain tensor was developed in terms of a Lagrangian curvilinear 
coordinate system which at all points remains parallel to the de- 
formed plate (see appendix A). The treatment of a developable 
plate problem in the Lagrangian coordinate system lends itself 
readily to inextensional bending theory, but becomes a major diffi- 
culty in dealing with a flat plate with extensional plate theory. 
Thus, the use of the plate strain tensor derived above assumes 
that the deflections are small enough that an Eulerian coordinate 
system can be used. Even though it is recognized that error is 
introduced by the assumption of an Eulerian coordinate system, the 
plate strain tensor does provide an improvement over the von Karman 
two dimensional approach as discussed above, and should lead to 
less error for any practical plate problem. Indeed, if point 
strains are sufficiently large to invalidate the above assumptions, 
it is probable that plastic deformation would also invalidate 
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the assumption of an elastic body and a non-isotropic or plastic 
theory would then become necessary. 

In reference [ 1 2 ] Bleich gives an excellent historical 
sketch of the evolution of the von Karman equations and subse- 
quent strain energy approaches to the large deflection plate 
problem by Timoshenko, Marguerre and Trefftz. In a series of 
papers by Levy [ 1 3 ] , [ 1 4 ] , [l5], [ 1 6 ] and [ 1 7 ] , solutions of a 
theoretically exact nature were given to the von Karman differen- 
tial equation ['ll]. In Levy's solutions, deflections and normal 
pressure were expressed in the form of Fourier series and sol- 
utions were obtained for various cases of loading and support. 
Although the method is involved and laborious, it has been recog- 
nized as one of the more accurate solution techniques available 
under the von Karman theory. However, the only plate solution 
known to the author using finite deformation theory incorporating 
the associated third-order elastic constants is a preliminary or 
exploratory study of a circular plate by Borg in reference _ 7 j. 
Although the solution in reference [7] is close to Timoshenko's 
solution of a similar problem contained in reference [ 1 8 ] , the 
approach to the two solutions was vastly different. Rather than 
using the classical approach of Timoshenko, Borg's solution is 
obtained from the non linear stress-strain relation developed 
in appendix B. This would tend to lend evidence to the validity 
of the Finite Deformation Theory approach to large deflection 
plate problems, particularly where both physical and geometrical 
non linearity are required. 
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Another non linear large deflection thin plate approach 
using the non linear inextensional bending theory is described 
by Borg in references [5] and [8], In this approach, a theory 
analogous to the Murnaghan non linear elasticity theory is 
developed. In the non linear thin plate inextensional bending 
theory, it is assumed, essentially, , that the large deflection 
is one which occurs without the development of membrane stresses 
and a large deflection form of the bending deformation relation 
is obtained. [8], The non linearity occurs in the second derivative 
of the deflection terms and arises as a consequence of the analo- 
gous behavior of a thin plate to that of an elastic body in gen- 
eral, and in a sense, bears the same relation to linear thin 
plate theory that the elastica (column) solution bears to the 
Euler column theory. 

If the deflection of a plate is not small, the assumption 
regarding the inextensibility of the middle surface of the plate 
holds only if the deflection surface is a developable surface, 
therefore, there is some question about the generality of the 
inextensional theory. Although the non linear moment-curvature 
tensor equations analogous to the non linear stress-strain rela- 
tions, as set down by Borg [8], provides a possible solution to 
the non linear plate problem, it would seem that the most general 
differential equations would result from the non linear stress- 
strain relation, where the classical equations of equilibrium 
incorporate non linear stresses, and the Compatibility equation 
results from the non linear stress-strain tensor equation. In 
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other words, the large deflection plate problem should be 
treated as one which is non linear both physically and geometri- 
cally, and incorporates the third order elastic constants. 

Borg [6] has shown that the order of magnitude of the three 
additional elastic constants resulting from Murnaghan's stress- 
strain relation are of the same order of magnitude as the elastic 
constants Or & E (see Table 1, page 5). As a result, their 
inclusion in the large deflection plate problem may be signi- 
ficant and would certainly warrant the derivation of more gen- 
eral differential equations. 



- 17 - 



DERIVATION OF PLATE COMPATIBILITY EQUATIONS 



1 . General Discussion 

Borg, Hoppe and Kopchinski derived non linear strain 
compatibility equations in reference [ 4 ] by: 

a) First, assuming = 0 

b) Second, neglecting entirely the requirements on 6 h 
and simply utilized the non linear relations for 6x f and 

as obtained from the non linear stress-strain relation. 

(6) T - Al, E 3 +2G -r^ + if- 2 mI*J)E 3 +2ml.yj +■ n Co^? 

The major limitation of the derivation in reference [4] is that 
it uses a stress-strain relation which neglects compressibility 
terms and other terms of the same order of magnitude as those 
retained. The more complete stress-strain relation, as derived 
by Murnaghan (see appendix B) is: 

(7) T=Al,ET 3 + 26n + [(/->;!, ' 2">ljE3 + + "coj! + 

In addition, the derivation given in reference [ 4 ] assumes = 0 
as discussed on pages 13 and 24 which in effect says that the 
deformation distribution across the thickness of the plate must 
be constant and ignores any Poisson's effect across the thick- 
ness. This is not an unreasonable assumption for very thin plate 
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but is a constraint that is open to some question for plates 
of finite thickness. (See the discussion of the order of mag- 
nitude of on page 24 following equation (19) ). 

Using an approach similar to that used in reference [4], 

Borg derived improved compatibility equations, the results of 
which are reported in reference [6], In the improved derivation 
Borg used the second form of the stress-strain relation that 
includes compressibility terms of the correct order of magnitude. 
Borg's derivation reported in reference [6] is verified in the 
derivation shown later in pages 31 to 37 • It should be 

noted, hoAv’ever, that Borg's improved derivation assumes r Q 

and thus ignores any Poisson's effect across the thickness of 
the plate. 

Intuitive reasoning A\rould indicate that although the 
large deflection plate problem is neither truly plane stress 
(as implied by the von Karman theory), nor plane strain, with 
the exception of the very long narrow laterally loaded plate, 
it more closely approaches a condition of plane stress than 
plane strain. In linearized theory, a plane stress assumption 
would require that Y xZ . = Yyi ~ ^ x 0 and " jTy (^6 X + e^) 0 

Thus it would appear that a more valid derivation in the general 
case, for plates of finite thickness, would be based on 0 

as shown in the plate strain tensor derived on pages 11 to 13. 
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2. Derivation based on generalized plate strain tensor 



with 




By setting = 0 and solving for in terms of 

oi 



w, x, y, u, and v, we can obtain a relationship which is in 
accordance with classical large deflection theory and which 
will be valid in the majority of cases (one exception might 
be relatively thick shell pressure vessels subject to very 
high external hydrostatic loads, i.e. very deep diving sub- 
marine pressure hulls). 

Thus, solving for the terms of the stress-strain 
tensor equation 



( 7 ) TsXI. E, + 26>? +[(< "A)!. 1- 3 - 2 (ni + X-G)Ii>J +• " 



+ n don +4-6 J 



we have (retaining terms up to the order €. ) 





J- <3* 






) 




( 



o 



0 0 



o 



o 



o 
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x. 

or, retaining terms to the order 6. , 



Cofactor matrix of the strain tensor 







^33 





Jn 



+ 



7 

I 3 ' 

)>' 



XI r /f3 



7 33 

^ 13 





7 33 




yj* 


?« 


^12. 


7' 3 




T- 












^13 



-h 
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Retaining terms to the order of magnitude € , we have: 

0 






w_ 



(10) 



Co 





± 


w ci w 


i&w \ 


r ■ 




o>X 


\9* 1 



(11 ) 



If 



l> 



*/" 7' 1 ^ 

j?U ^21 >/!3 

7 »i V 51 



> 7 »» fu 

fzi 

\ 7 ». *|» 1 



/ 





1 ^ 

** dj 

du/ 

< 9 X < 9 ^ / 

0 0 




lj and Xi contain terms of order higher than £ 
and are neglected. 

Thus, A^riting out the stress-strain tensor equation, we have: 



( 12 ) 



T = 



l<* 


T*\j 


*X*\ 


Tp 


6, 




[ r zx 


Bj 


w 



A 
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The equations for the elements of the stress tensor can now be 
written 



(13) 



(14) 



(15) 



(16) 



(17) 



( 18 ) 



6x * 





+ 2 & 










-h 




Z. 



% "a 2 

*f[(t iH%T] -f (fe)‘ * a 




d* 





+ 



'IL+2S.') 






cW ^vy 






G- ^ 
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Assuming 



0 as stated before, we have 




solving for &BL 

(19) _ ( A \{*J± 

V*+2G/\^X 



fa + m + 2 Gr) / 

2 (x+ 2 .&) [fa^^Uy 



It should be noted that — ^ — 

M2C- 



_*L 

/-v 



which for steel 



where y = .3 is about .43. 



Thus 






contains terms of the same order of magnitude as 



the terms retained in the basic derivation, ( 6 . ), and is of 

the same order of magnitude as the remaining terms in the first 
invariant of the strain tensor, and should, therefore, not be 
neglected . 

Substituting the expression for into the expression 

for and , we have 



(20a) 
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Rearranging terms, we have 



(20b) £ = 



+ 2 G- isi - XL ( * 2 l \ 1 4 2 G- ( $*L) 

** + T V s */ 



(21b) 



V 



^) ( * + ^ v )-[l(^) * f(a. 

+ ^t-; 



H-fJKfcAdf/J 



since 



>> 



?v -*-2& 



y 

i-y 



and 



A = 



y E 



C i+y)(i-2.y) 



2 6-y 

i -z y 



we 



can write the above expressions for <£* and in terms 

of the engineering constants 



<20o) 




+ ZG ( £u\ — H / &&) 

V <W -t l *<j/ 


(2.C 6 3 =(^rM^)-| 


m , 


'Z i 


+ 2 G /iv \ - J1 

V sy -4- 


(If) 1 



-2,^ V 1 






V a* 
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If we now define coefficients as 




(24) C = 2 6. 

(25) D = - fl. 

we can write 



(26) 

(27) 

(28) 



= 



"Zocij — 



a(I^) 





^ W ^vV 



Since we want to ultimately arrive at an equation in the form 



-£ < =^ 6 4 



which satisfies the stress 



function relation vV = ^ +• £ — rr" j. 



$2£ 



F 









where F is defined by 6 X = J 6i< - q ; r ■ 

J ’ J 3 x 1 1 J 3n3. 



q>* F 

o> Xt3*j 



as developed in the classical theory, we jjerform the indicated 
differentiation . 
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( 29 ) 



( 30 ) 



(31 ) 




Thus 



( 32 a) vV = + 0" - AV'(iSg) + B<f„ 






**/ 







Z ( D-c) 






<^w 9 w \ 

dx 



where da. = v‘ [($?)+(^f] 



1st invariant of the large 



deflection tensor. 

From the identity J-^ (§j) = ^ ( f* +- &) 



the tliird term goes to zero and we are left with 
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Carrying out the differentiation in the last two terms, we have 

C L Z (^) + Z % + Z (^) +2 ^^ **' *J l 2 <^ *3 MJ 

(" d^w _ / y w N* 1 
— ~2C v.<Sx^^/ _J 



, -_j /$V\ l , ^)w^w , /d l w\ Z . i!w , , Jwd’w /,: 

and 2 D[_(|^)+ 57 dx i + (^V f dij® Oy* f 7 T ^ Sx 1 ^ t 



^w y 

' V J 



If we consider the second invariant of the curvature tensor 
j _ Q) 1 W _ / ^) 2 yy N 2 - 



O ic 5>x 2 yo>x^(j / 

and considering the exjmnded form of the 1st invariant of the 

large deflection tensor ^ pu> = V = 2 [(0)V4&) *^7,)^ 



Ow ^ 3 W Q w «^ 3 w Ovy < x) 3 W 

£>X ^ X 3 c * >( ^.o?t^ 3 O/. °^ 2 ed j 



By adding and subtracting 



£ D 



/_<^w_y 

) 



we can write the last two terms as — 2c^c 4- D ^Iilp -t-2.0 ^ic 
Thus, the equation for V F can be written as: 
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(33) V“F = AV^rf?) +(B,°)d,u, + HD-O&i 



In order to get rid of the term A V + ^ we go back to 



the stress relations 



<26) • A(t B [( i?A(!^J + C £ * d(&)\ c (§?) 1 

< 27 ) M>( 4 £)‘ - c(|-:/ 



Adding we have 



w 






solving for + we have 

( 35 ) + = / 

O/. C^L. ( 






z a + cy 



4 A +■ 5'j 



- (*■& + P^o)[i~ 



Now taking the Laplacian of the above expression, we have 

(36) =(d;c)[ v ‘( 6 * t6 )) -( 2 BfD -cAt L » _ 

Taking the Laplacian V 2 {'£*to l ,') > by adding and subtracting 



„ ajjT- 



* 



we have ('g+-Z l'**» + **<$>) +( 1 ^ ' 

. ^a-3uj 3^ / V ^ A J.-A . 



< 3 q F , oOl^F , 3 F" 

^ V V 



;3 V 7 _ ? y*p y^ 



V 4 F 
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Thus we have 

(37) V‘ (I* v %) =(ii7 c ){v‘F - (2B *OA t .J 



Substituting this expression into 

(33) V'F = a[v‘(& 3-| p] +(6.D)J |L0 «(d-c)4 k 



gives 

(34) V V - 4* + 2(D-c)^ J 



since /\ = 2 (S- F 

1 -v> 



— 122 ( i - Z >0 

~ & O-vO 



c- 26 i D'-| 



we can write the above equation in terras of the elastic constants. 



Since 



( B+D) 

ZA ±S 
A +-C 



A (2S + Pi-c ) 
2A+C 

= ^+! 






-1 ; 



and 



E(P-C) =-Z(^+26r) = 



we have finally 



(35) 



V*F = [ 2 ( 1 -^) - £ -2.Crv\S, LD -(F+I )('4c-4Q)4, 
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3. Derivation of Plate Compatibility Equations based on the 
assumption ^ = o 



For the assumption ^ = O , we have the identical 
expression for the strain tensor as derived by Borg, Hoppe 
and Kopchinski [4]; 




i 

Z C)/; 



\ <? w 





Thus, the terms of the stress-strain tensor equation 

(37) T = 41, E 3 + 2 Gy + (/-A) l! - 2 m I 2 ] E 3 + 2(rr> + A - g) I, ^ + n C«f? 



are (retaining terms to the order e* ) : 



(38) i. * (£+!*) + }[(§?/ +(p‘l 

° 9) -*[(&/♦ (t?/] 
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( 40 ) 




(41 ) 





Terms 1,^ and I, contain terms of the order and can be 

neglected . 

Substituting the above expressions into the stress- 



strain relation gives 



(42) T =r 









1 \ x 



+ 2G- 



JL / i l 0 vu \ 

* <ix 






J- / W \ 

* K57) 




o 

\ 

o 



6 

0 

1 



!<• sy 4 S* + ^ ZjJ 






o 

/ 





O VO <^UI/ 



a 



cP W ly/ 

/^y \r 
l^>x/ 




O 
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The equations for the elements of the stress tensor (equation 
[42] ) are as follows : 



(43) .*{& + !; ![(!?/*■ ? 4 



, (^\ Z 
Six _ 



x /£*)*] - Hi-" 

* H**/ +{*•)/ J 4 



(44) ^ 



^ u , i j_ 



* *[% + df j‘] + z[(pjm%)i - f(ir) : 



<45) ^ + j [(W'<p 1 + 



(46) 



r *y = F - 



M/ 



(47) r„ . r, x - & 



i^W 

£>x 



(48) ? 



y 2 -- 



6- 






Since, regardless of assumptions, we want to eventually have an 



of a stress 


function 


F* or 


v 4 f 


^lr 


or 


o> Z 4* __ 


^ 4 F 








^F 
£ x"- 


or 


o) X 6u =. 





_ ^ F - 
- £xo>^ 



sr 



- £ = 



<S* <5 



* 
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we can carry out the indicated differentiation and add the 
resulting expressions 
Thus, we have 

iAn . / &< +$*) + 5 \<£ , ?r ^ (&K 

(49) + 

+ p r ^ . m £ -J1 A* 






’3 



<50 # - (t^) +\ & m ♦ **& m) 



^ (^-f + ^- — f (-) 1 - r] 

+ 2 £x 2 L\ d >^ 



(51) i 



<$* 



Tf v d( j 



n d* 


(— / 


4- ax 2 


Ux/ 




/ 4> w- ^y/ \ 


V 


Vc^x / 



Thus 

(52a) 



* 4 p - ^ - 2 fS; + 0 - * v \% + f^) + ( t=) v W/Hlf/] 

+ z& [£>(*) * fy(W 0$ *£)] 

+*& (tr -& (^;- ig^is $*)] 



_il Til \ x 4- £- 4-2— {£»&*) 1 

4- L l<3<^ J ^X* \^X/ ^X<^ 
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From the identity 




^ /du_\ _ 

\ e)x / 



and the definition of the 1st invariant of 
Tensor ^ + ^) 



the large deflection 




ve have (carrying out the remainder of the differentiation) 



(52b 




A±m 

Z. 



&LO ~ 




<^X ^ 






-I- I ^ w c) 3 w 




From the second invariant of the thin plate curvature tensor 

Jf ^ x w yw /«^*W 'J i 

c5ic — ax* (d^J 

/^W 

by adding and subtracting ve ^ ave 

(52c) V^F = AV l (f*-^) +. [W|4.u> -•‘tft^c -a (^, lb + 2^ jc ) 



Combining terms ve have 



(53) vV = XV*(fe+%) + U( - ^ n?) ~ - n] 4 ld -(^+-f)iac 



If we now asstime 6i 

>(S^tp= - 



= 0 , we have from equation (45) 

( 7(4 m + 2<5r) T (\>w \x , /^wnM 

a — LLdx/ J ^ 



( Vt- 01+2 6) ^ 



JLD 
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Substituting this expression into equation (53) above, we have 



(54) V F 

( 6* =■ o) 

(This form agrees with that derived by Borg [6] ) 

If we now make no assumption on and work entirely with the 

elements of the stress tensor 6* , 6y and , as was 

done in reference [4], we have, by adding equations (43) and (44) 

(55) 6 ** 6 ) + (* 4 )[(§*) V(f*/J 



Taking the Laplacian of each term, we have 

(56) V ( <5x. + 6^ ) - Z(~\+G) V +rri r ^‘ Cr ' ^ ^~/iud 

By adding and substracting o. ~ > we ■* iave 



(57) 






r-, z ( j-<^Y \ 

Thus, solving for v <3*. ) , we have 



V 

(58) \7 5 /~ - S 7 -— - (A + M+2& <A 

iO+G) 2(>-h6r) ‘ =/|LD 

Substituting equation (58) into equation (53), we have 

(59) V * F = ^ 74 Ft (*V**"**' w 



Combining terms and reducing, we have 



(60) W V _ -2(At6) f_ M **<**&■-%) 

V h A + 4& \ z (h +&, ) + -^T i Ju “ V 
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1 



or since 



_A_ _ _L- 

\-V 



~x 



2& y 

l -xy 



A 1 - 6 - = — 

1 - 7 .'/ 



2 




i-V 



y 



_j\ 

i (T ^ 



we have 

(61) vV 



I-/ 



+ i{ 






-(&)(** + ’$:) die. 




Reducing and collecting terms, we have 

(62) W = - £0-1') +- 

( - <3 , No restriction on 62 ) 

Equation (62) is identical to that derived by Borg [ 6 ]. 

Thus, we now have three new compatibility equations under the 
assumptions indicated that can be tailored to a particular 
problem. 

(35) V V -[?('-*■->-£ - 2 &*>]*„ -(yiX+fr*i) 4 . 

FOR USE IN PROBLEMS THAT APPROACH 




PLANE STRESS 



(54) VV 



" (j* + ^ ii_r» 



(46- +- ) ^ 3C 




FOR USE IN VERY THIN PLATES WHERE 
POISSON'S EFFECT ACROSS THE THICKNESS 
CAN BE IGNORED AND WHERE THE PROBLEM 
APPROACHES PLANE STRESS 



( 62 ) VV 



m 0 -2>0 _ £L 

2 . (i-O 4 




0-0 + 0 ^ lc 



= 0 

^ 2 . 

No restriction on 




FOR USE IN PROBLEMS THAT APPROACH 
PLANE STRAIN OR WHERE 6z MAY NOT 
BE NEGLIGIBLE AND = 0 IS JUSTIFIED 



And finally, for comparison, 



the classical von Karman equation 



(63) 



V^F = -£■ 4xc 



- ( y+i) 2 & d ic . 
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DERIVATION OF EQUILIBRIUM EQUATIONS 



In the classical linear theory of elasticity, the equations 
of equilibrium are 





div T" 




= 0 




or 






(64) 


^ix _ 4 _ ^ 




4- /=>*. Fk = O 


(65) 




t ^ 


+ /°* F i = ° 


(66) 


o>y. u 


<2>z. 


+ Fi = & 



However, for the large deformation case, one must remember that 
the above equations are referred to the original coordinate system 
and that the curvilinear coordinates a, b, c, are identified with 
x, y, and z so that when the forces are projected, the changes in 
position of the points of the body due to its deformation are 
neglected. Stated another way, in the linear theory, no distinc- 
tion i‘s made between undeformed and deformed values of the magni- 
tudes and positions of the elemental areas on which the stresses 
act. In other words, in projecting the forces, the rotation which 
an element of volume experiences as a result of deformation is 
neglected. This assumption is open to question and is far from 
admissible in the general case. In the general case, it is nece- 
ssary to take into account the fact that differentiation is with 
respect to a, b, c (the coordinates of the points before deforma- 
tion). Since the nature of the simplified linear equations of 
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equilibrium rests on the assumption of small elongations, shears 
and angles of rotation, and in view of the interaction of these 
quantities in the higher order equations, whether or not the 
non-linear terms can be neglected depends not only on the mag- 
nitude of the terms of the strain tensor, but also on the com- 
paritive magnitude of the corresponding terms of the stress ten- 
sor. It follows that the smallness of the angles of rotation in 
comparison to unity is not a sufficient condition for lineariza- 
tion of the equations of equilibrium. It is also essential to 
know whether the stresses which are multiplied by rotations are 
large in comparison with those stresses which enter linearly into 
the equations. The problems of elastic stability or of thin plates 
with large deformation are cases in point. 

The curvilinear form of the equations of equilibrium 



(67) 


^ 6* L 


i ^ (- b*. 


q_ fea. 


4- = o 


<£> a. 




c. 


(68) 


drab 






+ /°* F b - ^ 


0 <x 


$ b 




(69) 


dr* c 




d 6c 

^ c 


-4- /Oc F c - d) 



can be transformed to the cartesian coordinates of the points 
of the body before its deformation. Novozhilov [2] has derived 
a vector form of the generalized equations of equilibrium as: 
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where £n (> and 6n } are vectors, 



= 

5* 






•>V 



0 



;£h* = 

6z 



Ratios of area of 
elemental sides before 
deformation to area of 
sides after deformation 



J 



V 

V 



ratio of volume of element after deformation 
to volume of element before deformation 



= D 





i ‘ 


Ou 






*3 


3>Z 


D - 


Ov 


/+• 

* 1 r 






d>/. 


$ w 
** 


| i 


and 


I 1’ *2 


and 1^ = 


1 St, 



= | J I = V I + 21 , SI, 



Xu Z 



Resolving the vectors 6n t 6n 1 , 6>n 3 in the directions i, .j and k, 
we can write Novozhilov's equation 

< 71 > k [4f ( ** +T ii r f r ? £ ^)/ 

4 5i [jsz + + 



y R = <3 
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( 



Projecting this vector relation on the x, y, z axes we have: 



( 72 ) e,, - Ul ) e , t +W) ) v * ] 

ECl * + (x^*£ + J 

+ Ij[0+ex)r 2 ^(i€x r ^ 2 )^*+(ie xt+ ^) 6z] + "V F * = 0 

(73) ^|]i e x t 4^ 2 )^ x * + (i+<?>) r x f +(je^-£j*)r K t J 

+ 5 j r (i. ex ^ + +-( 1 + 6 ^) 6ij, j ^ 

+ + £J «) Tjw t ( i+ep ] + V Fb ~° 

(74) §- % [(ie Xi -^)6? +(ie^+«*)r£ +(i+ejz£] 

+ ^ 3 x ^7 +■( l+ 6 z) 

+ & e « 5 * +(i e^ a -hcox) + ( i+€i) 6z] + V~ Fc = 0 



where 



e« = ^ 
Cx 


i e i 


_ 




p _ sLw. 4. Sv 

) e ** ' ^ + 57 


0 

) 


«« ‘fS+S* ; 


e i* - & *• 


4 

) 


-2 ^ 

57 




*«»* 51 


- 

»?x ) 




» 

) 


. * si 6x 
<^X * <5j« 1 l-ETx 




r * _ ' 5 * 
**» ' «5 a 


Jh- . 


r v , -*. r-„ 

Kl S* 1 +£i 




/* - *£. M- . rl 




» + tr<j, ) 


) 


'* E ^ 1 f* 


^ ,+£i 


K Su 


?V 

I rex > 


* S? 69 
6 i ‘ ^2 /*•£> J 


~* a , Tit * 

CiK /re x > 


i r*- 
■ 3i 



fT x , £<j a*«l = Elongations in fibers in directions x,y and z 
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The above values of normal and shear stresses are not, strictly 
speaking, stresses. They can be called stresses referred to the 
dimensions of an element of volume before, not after, the deforma- 
tion. [2] 



Writing the non linear form of the general equation of 
equilibrium in terms of u, v and w (the displacements) we 
have : 




(77 ) dw/ f h 












c)x^ 



ii/v <^x 
* 



h 4i/ * <>x' 6 ' 



A *" 



+ 2 v v. 



- * 



* * Ss 4 






+ -2 £>u.di ^ 



>? 

* \/+ _* 



dV / ' 

5p*J 



* A* 

o)* 



+ = o 
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Txie above relations clearly show the interdependence between 
the rotation and the stresses in the equations of equilibrium 
in the general case. However, the above equations are in gen- 
eral far too complicated to be of significant practical value. 

Since the basic assumption of this thesis is that the 
large deflection plate problem is both physically and geomet- 
rically non linear, the classical equations of equilibrium can't 
be used without a further analysis since their derivation is 
based upon 

(1 ) The assumption of plane stress 

(2) Geometric non linearity without regard to the physical 
non linearity. 

The classical equilibrium equations (appendix C) are valid 
approximations of the more general equations (75), (76) and (77) 

and thus, the geometrical portion of their derivation remains 
valid for the purposes of this thesis. Thus, the equation 



(78) M x —9 — 
ex’- 







remains 


valid, 


and the moments 


still are given 


V 


(79) 


M x = 


J <Sx i? c/h. 






A 


(80) 




f % 


M - 


-J 


(81 ) 


** 


- [ T**zcht 




& 






per unit length and twisting moments 
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At this point, however, a departure from the classical 
approach becomes necessary. Instead of using the linearized 
stress relations to obtain a relationship between moments and 



deflection w, we will use the non linear form developed in 
appendix B, i.e. 



(7) T= *I,E 3 + 2G>^ +[(*-*) I,*- 2m I*]e 3 + + nCoij +-4 & y ^ 

Using the analogy between the stress-strain tensor equation 
and the moment-curvature tensor equation as suggested by Borg [5j 
we have: 

(82) M = ^l c E 3 +2Gt%+[U-X)d,c-2m4^E 3 l-2.(rn+?i-G)d K $ + n c 0 tf 4-f^ 



M - 

and 0 ^ 
(83) 



is the 


moment per 


unit 


/ 


■ M a i„ 




Mba. 


Mbb 


-(2b^ 


\-GaC 


-ObC 


z < 


is the 


curvature 


tenso 



<f = 



/-L 



rca. 










/tffJ 

\ 

i^V/dV 

5? 5Sv> afe 1 JoAc 



/ o>w \ 
dcUc/ 



| c)\v J t w . 

d4 l ^o^c 

tJbic dc v 
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If we now use the classical large deflection theory assumption 



that the 


slope 


is small 


and the 


square 


is negl 


igible 


to unity, 


ana 


assuming 


Eulerian 


Coordinates, we have: 




i-i 


1 

s 


l \ 




/& 




c>V \ 
\ 


(84) 


_L 

V 


f 


V 


- 












r *J 


1 

V 








/ 

aw 



At this point, it is necessary to consider what actually happens 
to the plate when we make the assumption ^ O . This assump- 
tion, in effect, assumes a deformation gradient through the thick- 
ness of the plate (Poisson's effect) which alters the limits of 
integration in the expression for the moments per unit length. 
Thus, to be exact, for the assumption o we should write: 
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where £ is the change of the neutral surface due to the deform- 
ation gradient ^ across the thickness of the plate as shown 
in figure 2 below. 
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Deformed Locus of Centroids 
In the Middle Surface 




Neutral Surface 




5 = 



4r+Jt 

Z. 



From figure 2, it can be seen that the total thickness is 

essentially not altered by either the assumption ^ =£ o or^5 !—q , 

o>H ' dZ 

but that the plate deflection W is altered in the case of = o 
by an amount £ > which in effect, requires an additional constrain- 
ing force F(£) such that the total deflection is now w» VV(x,tj)+w(<jl)* 

Thus S is in effect, a measure of the error of the assumption 

= o 

A complete investigation of the effect of the S factor 
is beyond the scope of this thesis, however, it can be seen from 

equation (84) for the curvature tensor, that the curvature elements 

d W 4- 

containing •— — - are all of the order of magnitude € as dis- 

cussed on page 13. Accordingly, to the degree of accuracy of 
the other equations of this thesis, we can write the curvature 

tensor as _ 

/ jv 1 ... r\ J V y^ 
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From equation (88), we can see that we are left with the express- 
ions for curvature that are identical to those contained in the 
classical theory in two dimensional form, where the third dimen- 
sion terms are assumed to be zero. 

From the curvature tensor [equation (88)] we have the 
following relations 



( 89 ) 4,c 

(9°) 4 



/_! J , J_ \ - 

w* ^x 1 - y *■ 

4“ ~r ~ (t~) +■ tcrmi involving 4 “ , 7" An d -L_ 

T * J h r *z r a z 

<^\y cl^W _ / 
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can 



contains terms 
be neglected 



of higher order than desired and thus 
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= Terms of higher order than desired ( € ) 
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at 



— terras of the order 6. and can be neglected 






terms of the order 6 and can be neglected 



Thus, we have, retaining terms to the order 



K ” 



$I*y o>Jt/ 



It should be noted at this point that retaining only those terms 

a. 

to the order £ reduces the Moment-Curvature tensor to a form 
that is analogous to the linearized Hooke’s law T* >I,t 3 +26>? . 

Thus, we can see that the classical assumptions in effect force 
the equilibrium equations into a case of geometric non linearity 
and physical linearity where the third order elastic constants do 
not appear. This fact is born out by the fact that the third 
order constants 1, ra and n are multipliers of non linear terms 
and cannot appear under the classical assumptions for the equi- 
librium equation. 

Inasmuch as we are now left with a linearized stress-strain 
relation and its analogous Moment-Curvature relation, the classical 
equilibrium equation derived on the basis of plane-stress becomes 
valid for cases where the actual case approximates plane-stress. 
Thus we have, from appendix C and reference [18] (neglecting body 
forces ) 



(94) v’w = y + 

(PLANE STRESS) 



t_/ 

P { Ju' 
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Proceeding with the development for the assumption = o 






we have, from e= -y- in the two dimensional case and the 
analogy between tensor equations 

6% K 



(95) 

Thus 

(96) 

(97) 

(98) 
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Substituting equations (96), (97) and (98) into the expressions 

for moments per unit length, we have 

t/z 
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Performing the above integration, we now have: 

<102, M x = 



(103) M 



- i! I ^ + ^V ), i 

if - TzLM*x' + a^ 26 ^J 
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(,04) + ^ 



Differentiating equations (102), ( 1 03 ) and (104) to obtain 



£ x Mx _ 2 + 

<2> X* 



, we have 
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adding, we have 



(108, _ 2 ^ + ^ = 



e x r 



_±! (*+;>&) V 4 W 

/2- v x 



Using the relationship Th-2-O = ££ £ 4- 2 6- = 2£_(±lE) 

i-Z>> I-2.V 



and G - fc 



Z(\+9) 



f ve have /\ ±2& - *=• £<->0 

I - 2.v0 



Since the flexural rigidity is D =. ILt — - 

12 

We have finally 

( 109 ) - 2 — ~ J> ilzll V 4 W 

Sx 1 Xw4i/ ^u 1 - fl-2/) 






Q-20 



- 51 - 



>* 



The equation of equilibrium for the assumption 
now becomes (neglecting body forces) 



S|w = O 
<=*z. 



( 110 ) 
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or from the stress function F defined as 



(Hi) 

( 112 ) 
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We have 
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As a check on the above equation, since the basis was a linear- 
ized tensor equation, we should be able to obtain the identical 
expression using the linearized plane strain equations resulting 
from the linearized Hooke's law 

T = ^\I,E 3 -*- 2 G- where for plane strain 

we have the elements £ ^ =- Y« Z - Vx! «* = o 

d> £ 

Thus, we have the expressions 

(H5) 6 x - "A(e* 4- Z& ex 

(116) 4- £j 0 4-26-e^ 

(117) r x ^ =■ T r - 
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Since, in the linear relations 






2 Z. 



where 



a = - 2 .& 0 . 

0 * 
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= A = - 2 
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we have 


the expressions 
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Equations ( 1 1 8 ) , (119) and (120) are identical to equations (96), 

(97) and (98), thus we could have derived the equilibrium equation 
(114) directly by the plane strain linear relations. The point to 
be stressed in the derivation of equation (114) is that the classi- 
cal large deflection assumptions, when coupled with the assumption 
^ W 

— O leads to a plane strain equilibrium equation rather than 
a plane stress equation as derived by von Karman. 

In summary, within the limitations of the assumptions of 
classical large deflection plate problems, we now have two equilibrium 
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equations, the von Karman equation where the derivation is based 
on the linearized Hooke's law with strain a function of stress 
and the equation derived in this thesis with stress a function of 
strain. Thus we have an upper and lower bound on the assumption 
of plane stress or plane strain as follows: 



PLANE STRESS 



(Neglecting Body Forces) 



(94) V<W 



f, t_ /^F ? a a W \ 

D D l ~ 4 / 



PLANE STRAIN 



(Neglecting Body Forces) 



(114) V*W= + t 4- 

As previously discussed, the actual case will be somewhere between 

the limits of the above equations. 

I — 2P 

Applying the correction - — as a linear correction to 
the flexural rigidity in the data reported in references [14] 
and [l9] and comparing the results with experimental results 
shown in reference [ 1 9 ] we can see by figure 3 that the correction 
appears to be about the correct order of magnitude and in the 
correct direction. Thus, for this case, it would appear that the 
actual case is closer to plane strain than to plane stress. On 
the other hand, one must be cautious about drawing hasty conclu- 
sions based upon figure 3 since there was insufficient data available 
in reference [ 1 9 ] to accurately evaluate the boundary conditions 
of the experimental results. Further discussion of the implications 
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of the correction of Levy's solution and the divergence of both 
theoretical solutions from the Galcit experimental results is 
contained on page 62 . 
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1.6 




CONCLUSIONS 



In addition to the classical large deflection thin plate 
equations derived by von Karman, we now have three new compat- 
ibility equations and one new equilibrium equation based on the 
assumptions indicated below: 

A. COMPATIBILITY EQUATIONS 

Classical von Karman Compatibility Equation [ 1 1 ] 

(a) v‘f = -E4 = -(y+i )2 g &^ 

Based on = 0, linear stress-strain relation and plane 

stress 



Borg Compatibility Equations [ 6 ] 

(b) V*F = -(.(,*%) -(44t2)^ t 

Based on 6z. = 0, = 0 and a non linear stress-strain 

relation 



(c) = 

Based oii^ 

es 

stress-strain 



f m ( 
L z 



~ 4 t> 4 j)^ic 



= 0, no restriction on 6% and a non linear 



relation. 



Compatibility Equation Developed in this Thesis 

(cO v'f = [ f ('-*>!) -(y+0(4&+ 

Based on = 0, no restriction on o>w and a non linear 



o 2. 



stress-strain relation. 
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B. EQUILIBRIUM EQUATIONS 



Classical 

(e) V W 

Based on 
strain is 



von Karman Equilibrium Equation [ 1 1 ] 

f jb O x W 

“ i> + X> dx 1 9x x £ 3 *- 

a linearized plane stress assumption 
a function of stress). 



) 

(Assumes 



Equilibrium Equation Developed in this Thesis 

, , «<• _ l-2>> f t fl-lf) _ p \ 

^ ; v w - D + p J 

Based on non-linear development approaching plane strain, 
but linearized by use of the classical large deflection 
assumptions. (Assumes atress is a function of strain). 



From the above four compatibility equations and two equilibrium 
equations, we can now tailor the equation to suit a particular 
large deflection plate problem dependent upon the degree of 
accuracy required and the degree and type of deflection or 
loading involved. In other words, the equations can be tailored 
as follows: 

(1) Equation (a), or the classical compatibility 
equation can be used where the problem is linear physically, 
but non linear geometrically (moderately large deflections), 
or where the degree of accuracy is not critical. 

(2) Equation (b) can be used for ver) thin plates where 
the problem is non linear both physically and geometrically, 
but where the Poisson's effect across the thickness can be 
ignored, and the assumption 6 b = 0 is valid. 

(3) Equation (c) can be used for thin plates where the 
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problem is non-linear both physically and geometrically, but 
where the Poisson's effect across the thickness can be ignored, 
and where the normal {62 ) stresses can no longer be neglected. 

(4) Equat ion (d) can be used for moderately thick plate 
problems non linear both physically and geometrically, and 
where the Poisson's effect across the thickness of the plate 
must be included, but where the assumption 6t - 0 still remains 
valid (Thicker than Thin Plate Theory). 

(5) Equation (e), or the classical large deflection 
equilibrium equation, can be used in cases where the boundary 
conditions and loading are such that the actual case is closer 
to plane stress than plane strain. 

(6) Equilibrium equation (f) can be used for problems 

where a plane stress assumption is not valid and where the 
Poisson's effect across the thickness of the plate can be 
ignored (~ = 0 ). 

Reflecting on the above equations, the question arises: 
why do the von Karman equations predict results that often closely 
approximate experimental results, but occasionally are signifi- 
cantly at variance? Particularly interesting is the fact that 
for the preliminary solutions Borg has achieved with equation 

(b) for circular plates with circular symmetry, the results are 

/ 

very close to those obtained with the classical von Karman equations 
(a) and (e) L T ] . In view of the apparent wide variation in equa- 
tions (b) and (a), one must ask - why? 
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If we neglect the terra in equation (d) involving the 

large deflection term 4\ld > we can note the similarity to the 

/ 

classical von Karman equation i.e. 



(a) V«F = 


-O-M) ICr 4xc 


(d) V m F - 

modified 


- (y+»)(^ Cx + 



From equation (a) and (d ^ m odif ied’ we Can See ^ 

third order elastic constant n is equal to -4G, the von Karman 
equation is identical to equation (d) if the terra 
can be neglected. 

As shown in Table I, Borg obtained values for n (corres- 
ponding to Young's Modulus of 30 x 10^ ) of -34.2 x 10^ . 

Also, Smith obtained values for Austenitic Steel of -40.0 x 10^ 
and for 0.6 Carbon Steel of -67 x 10 . If we assume an average 

of the above values of n = -47 x 10^ and the commonly accepted 
value of G = 11.5 x 1 0 6 , or 4G = 46 x 10^ , we can see that 
the actual numerical value of n probably does in fact lie 
very close to - 4G. Hence, we can approximate equation (d) as: 



V"F 2 [(§ +&)(t-2>>)]d,a, - (V*l) c 

J / 

The above equation still does not explain why the von Karman 
equation (a) provides valid results in the majority of the cases 
since we have only explained the right term and not the term 
involving 4 iu> . 

If we analyze the function 4iu> ~ + 51J* 3 ((.§* ) J 

by vector algebra, we can see that grad w = = V VV 

and hence grad w - • grad w =(^) +{^p = I Vwf 

but the maximum directional derivative dw 

d S m** 



= I 9 



rad vvj 
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hence, the vector meaning of o/u> is that it refers to the 
Laplacian of the square of the maximum directional derivative, 



i . e . 



V 1 




x 



But in the classical large deflection assumptions, the square 
of the slope was considered small compared to unity, therefore 
the Laplacian of the square of the directional derivative will 
be much smaller than unity. Thus, within the limits of the 
classical assumptions for large deflection plate problems, the 
function involving <^iu> will be negligible unless the slope, and 
hence, the deflection is very large. Stated another way, when 
the deflections become such that the term involving A ild cannot 
be ignored, the classical assumption, in which the square of the 
slope is much smaller than one, is questionable. 

If we further consider the elastic constants associated 
with diLD i.e. ~ ~ - XQrV ] an<J assume that n = -4G , 

vv = - (V+!) ^ 4x, 



we have 



From Table I, for Borg's and Smith's values of m, we can assume 
an average value for steel of approximately ra - -60 x 10 , or 

about -5.2G. With this assumption, we can approximate the com- 
patibility equation (d) as **’ ^1* Gr)(\~ 2.?) ^ ILD ~~(y A-G' $ Xc. 
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or, for y 



6.3 



V * F 



— .64 G^ ilo 



or 



* (o Gr 



V 4 F r -&f.64^y) +^.6^c] 



Hence, ve can see that the elastic constant contribution to 
the terras of the classical equation is about 4 times as great 
as the contribution to the term involving . 

From Figure 3, it can be seen that experimental values 

for the simply supported square plate tend to diverge from 

We, 

theoretical values above of about 1.5 for the data presented. 

Thus, from the hypothesis regarding the contribution of the term 
involving <^iu> and the results shown in figure 3, we can further 
hypothesize that for very large deflections ( greater than 
about 1.5) the terra involving ^lt.D begins to become significant 
and the von Karman equations lead to increasing error with in- 
creasing deflection. Below the divergence area, the ^, UI > term 
is not significant and the von Karman equation is very close to 
the more exact equation (d). 

The classical equation of equilibrium (e) is identical 
in form to the equilibrium equation derived in this thesis (f), 

with the exception of the coefficient involving Poisson's ratio 

O - 

. For a Poisson's ratio of .3 the coefficient is 0.816 
and for V = .25 , the coefficient is 0.89. However, as the 
deflections become very large arid the deformations approach the 
plastic range, we have a different situation. In classical 
theory, when an isotropic material reaches plastic deformation, 
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Poisson's ratio is taken to be \ [23], However the coefficient 

I - 1*’ 

goes to zero at 9 - £ and hence, equation (f) becomes 
Vw = 0. Reference [25] states that some investigators have 
obtained values of Poisson's ratio in the plastic range much 

greater than ^ (as high as 0.8) which leads to further questions 

4 - 

since the equilibrium equation (f) becomes V W = -f (F,w,x, y) • 

As a result, we are led to the tentative conclusion that equation 
(f) is only valid in the elastic range where the range of 9 is 
0.25 to 0.35, or that for deflections approaching the plastic 
range, we can no longer ignore the Poisson's effect across the 
thickness of the plate (the assumption =0 is not valid) 

and a three dimensional development similar to the development 

of equation (f) with 0 will be required. It should be 

d 2. 

noted that the three dimensional development of an equilibrium 

equation where 4 1 0 , as proposed, would be based on the 

d 2. 

assumption that stress is a function of strain as opposed to the 
classical assumption of strain is a function of stress. 

One final point to be noted in the discussion of this 
thesis is that in equation (d) 

<d) 7*F = - C 



the derivation was based upon the classical large deflection 
plate problem assumptions, but including terms applicable to 
plates of finite thickness, and with finite deflection theory. 

As a result, the development involves a compromise between thick 
and thin plate theory. Hence, we have a third category of plate 
theory which may be called "Thicker than Thin Plate Theory". 
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RECOMMENDATIONS 



The work done so far in the application of finite 
deformation theory to structures has hardly scratched the 
surface of a wide and virgin field of science. To attempt 
a discussion of areas requiring further investigation would 
be the subject of a thesis in itself. 

A few of the more obvious areas of interest to the 
structural engineer which require further investigation 
are : 

(1) To obtain more reliable values of the third order 
elastic constants 1, m, and n. 

(2) To obtain solutions to large deflection plate 
problems by use of the equations developed in this thesis, 
and for wnich there are experimental results available by 
which one can obtain an idea of the significance of the third 
order constants, and to determine the validity of the hypothesis 
made in the conclusions of this thesis. 

(3) Determine values of S as applied to the limits 

of integration of the moment per unit length equations in order 
to determine the magnitude of the error involved in neglecting 
Poisson’s effect across the thickness of the plate. 

(4) Through dimensional analysis, attempt to determine 
a physical relationship for the third order elastic constants 
as has been done with Lame's constants in linear theory. In 
other words, attempt to obtain a relationship between 1, m, and 
n, and the defining parameters of a material that determine its 
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thermodynamic solid state point (Temperature, Pressure, Specific 
Volume or Density, Conductivity, etc.). The long term impli- 
cation of the above analysis is the development of a combination 
Solid State Physics/Thermodynamics approach in the determination 
of material properties. A logical extension would be to include 
time dependent functions for purposes of analysis of creep, fatigue, 
etc . 

(5) Attempt to determine equations for critical stresses 
for buckling or instability using the plate strain tensor derived 
in this thesis with a linear stress-strain relation in addition 
to buckling equations using the compatibility equations of this 
thesis which include the effect of the third order elastic con- 
stants of the non linear stress-strain relation. 

(6) Derive a Compatibility Equation under which no re- 
striction is made on either 6i \ or ( 6z ^ 0 and-^ £ 0 ) 

for use in "Thicker than Thin Plate Theory" where must not 

be neglected (as in the case of very deep-diving submarines). 
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APPENDIX A 



Strain Tens 



or ) ^ 



In tensor notation, the state of strain of a body, initially 
unstressed where point "a", having coordinates (a,b,c) , then 

strained to point "x" having coordinates (x,y,z) } is given by: [5] 



(A— 1 ) 



STATE OP STRAIN = 



dx*dx - da* da 



Since dx = J da, where J* = ^ 3 * — - Jacobian matrix of (x,a) 

dO> b, O 



or 



_ / 1a ^2L lit 

\J^ <^b <$c 



da. di> dc 

AZ dE ^2 / 
^ do. dto £«-/ 



= a tensor 



and 



d x = d a* J~ 



T 

where vj is the transpose of the Jacobian 




^■oC da a O. » 



^ <i z ) =\£ 12 H 

M 12 -H/ 



and is also a tensor, we can say: 
( A-2 ) STATE OP STRAIN 



- da*(j*J- E 3 )J« 



z = c + vV 



where E^ is the 3x3 unit matrix 
Since X = CL -f U. > = b + V 

where X - X ( Ck, b, c) , (j = tj ( a, b, c) , and EsZ(ft, t ; c) 



A-1 



we have ^ 



+ iid ij4 

4 ) 0 . , db 



|+I5 



a Z 

<5 c. 



t^W 

lie. 



H 



ae 



dii <$* - 

5 x 5To 5" b , ^ oc 

du. ^ i = Q>u 

d4 ti C 



V_ ^ 
ciCL 



etc. (using chain rule for other terras) 
we can express Jacobian as: 



Tt> 



Q>U 

*>c 



(A-3) J = 1 g 



a* 



+ 56 



do 



Ovv 

<i<\ 



dw 

£ b 






Sc 



/ 



and 

(A-4) 



J*= 






1+ 

^OL 




*5a 


& 


|p^ 


^A/ 


sb 


<£b 


&b 






l+|S! 


S c. 


&C- 


£e 



\ 



Expanding j*J we have : 






( A— 5 ) J*J- 



[ 1 +2^ +i!:: r^+.^ T ^^+^^+l*'ij) 

L **■ WH^VJ [5b f *L + 3i<*b a*;>v> a*db/ [* j*.*c 






f Aa^u ,iwJw~) . iv iv 1 

dc. do. Sc5c io-^c 4 o.<5c.( /^c ib ^b*«- Sb 






a symmetric tensor 
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Writing above expression in terms of orders of tensors, we have: 

J - E 3 +■ Z ^ + 2 ^ = E 3 4- 2 (?,+>?*) 

where )^ ( is a 1st order tensor 

and is a 2nd order tensor 

^ is the combination 1st and 2nd order strain matrix 



and rearranging, we have 

y**. ^ac 

(A-7) 







J*J - ET. 



From the expansion of and the relation 7 - i (J*J-E,) 



ve have 



<5 a. 



-i- \ \ 

^ ' db <ia.y A wc, J 



( A-8 ) 



7 



L($a ±*v\ 

i Ub 

X /*<£ 

2.{* c <9«k/ 






\ 



±-(&X 

lUc r dW 



I /$v , <$*/ 

2. Uc 

dW_ 

dc 



+ Z. 



/[/ \ x i$y \ v . /^w\M r^ay^u . 4y .. dm/ j>v~l f Ou 4 m . j- <d>v <Jw~7 

[|UJ VJ ‘V/J [**• *•“ * fcib J ^ <j<xszj 



The above expression expresses the strain tensor in the Lagrangian 
coordinate system where the movement of each particle in the body 
is followed (i.e. the initial coordinates of a particle were (a,b,c) 
then the final coordinates of the same particle in the deformed 
position at (x, y, z) is [ x(a,b,c), y(a,b,c) z(a,b,c)] ). [5] 

If the assumption is made that the deformation is sufficiently 
large that the 2nd order strain tensor must be included, but is 
sufficiently small that "a" is essentially the same as "x", the 
expression for strain can then be expressed in the Eulerian coordi- 



nate system (i.e. conditions are expressed at each point in the 
deformed body, and a particular point rather than a particular 
particle is considered and quantities considered are functions of 
x, y and z ) . 

Hence, we have the approximation 




(using American notation for elongation and shear strain). 
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APPENDIX B 



Derivation of Finite Deformation Stress-Strain Relation 
I. The fundamental relation of elasticity theory connecting 
stress and strain is: 



(B-l) T = j 



where T 


is the 




1 6, 


7- 


T = 


t y 












pc*. 


is the 



J Sf J* 

stress tensor 

\ 

) 

6i ) 

compression ratio 



J = Jacobian 



det T 



Li] 



o) (Cl, b, C ) 



0 = $(Y() = energy of deformation / unit 

initial vol. = 4^ 

strain matrix = 



(B-2) 



7 

f = £ 

(see Appendix A) 



II. If the deformable medium is assumed to be isentropic, it 

is elastically insensitive to every rotation of the initial 
Cartesian reference frame and hence = <j> ( R* >; r) 

where R = the rotation matrix. 

Further, a deformable medium is elastically isentropic if and 



only if is a function of the three invariants 1^ , I, 

I 3 of . 

B-1 



and 



In reference [ 1 ] it is shown that: 



( B-3 ) 
(B-4) 
(B-5) 



11 . = e 3 = 

^ -r 



/' 

o 

I 0 



0 

1 

0 






= Co 



7 



r 



It follows, that for an isotropic medium 



= ft t| * M ft ^ 









or 



I f =|^ E 3 



7 

4 - ^ Co 

<^r. 



7 



III. Assume the energy of deformation per unit initial volume 

is developed as a power series in terms of the three strain 
invariants 1^ , I 0 and 1^. If ^ is assumed infinitesimal, 



1^ is of the order of €. 
1^ of the order fc 3 . 



, I 0 of the order fe 1 , and 
We can then say: 



(b-7 ) 9>(y) * = A. a,IafA 1 I, 1 + jlA,!,-!-"' 

+ 0 , 1 ,. + <-•••* \ 

+ C, Xj, +■ +••• + D, 1,1^ 

3 

Neglecting terms of higher order than C we have 



(e-si #(„) = a. +a,i, t iA, i , 1 fX. a 3 i , 5 + e,r 2 + 0,1,1* +j),rj 
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i 



+ .. I 



or in terras of & - +■ <P, + <Px + <Pi 

where <p a consists of terras independent of W 
<f> x terms linear in elements of ^ 

$>X terras quadratic in elements of ^ 

- * . 

and hence a linear combination of 1, and 1 2 , and 
a linear combination of I| 3 j lilt and X 3 

(B-9) 4 - All, v V * fi . Ii) + % I.U tD,r s j 



<Pi 4>x 

Differentiating with respect to ^ , we have: 



P; 



(B-1 0) 

*r 




Since 



(B-1 1 ) 



Let 
(B-1 2) 
(B-1 3 ) 
( B-1 4 ) 
(B-1 5) 
( B-1 6 ) 







+ x>. 


X 3 

*7 




iL = IT* 




/ 

) 


|h = Co> 

^ t 


/ 


- A £ 

f D, 


+ Ax I, tr 3 +- s, (r (E - 3 

Co Y7 


->?) * 


I,X f|. 




(At + 8 .) - 


- vf 

* ~ 0 +*Xt-AV 


Lame 1 


s Constants 


U&/< ) 


8, = ~ l/*- ~ - 2. & 


( y 


= Poisson's 


ratio ) 


Aj+C, _ 
2. 


i 








c, - 


— 2. vn 


l,m, 


and n = 3rd 


order constant 


D» 


n 
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Substituting the elastic constants into we have; 



(B-1 7) = A,E 3 + ^I|I=3 p i- n cop 

which is the 3rd order approximation obtained by Murnaghan [l]. 
The higher order relation between stress and strain 

* vsy; 

•* 



/Opt "T - / \ ”T ^ 

resulting from the generalized relationship T = -rr J I v“ 1 *■' 



is then: 

( B— 1 8 ) T- -<■>!, Cj+^6-p +(/Ii*‘-imrjH'j4--2mXip J J~ 



If the elastic medium is initially unstrained, 



= 1 and 



X 



is a zero or null matrix, therefore J* is a rotation matrix only. 
If T is the initial stress in the unstrained condition, the 
above tensor equation becomes 

(B-1 9 ) T 0 -* A. E 3 

where T q is due to a hydrostatic pressure or tension. 

For an isotropic medium, we can say = — p (i.e. 

initial hydrostatic pressure) or T q = -p Q E^ . With an 
initial hydrostatic pressure of -p , equation (B-1 8) becomes: 

1 B-20) T= $ Crj-f^lTj +Al,f 3 + 2<»)7 4- (A I, -2 i^I^) E 3 4- + n coy^] * 



For an isotropic medium, the coordinates of the stress tensor 
are furnished in the final rectangular Cartesian reference 
frame in which the coordinates of are furnished by the 

elements of M = by the elements of the matrix /f* M 



where J r 1^ M and J~ * J" = M = 2 ^ 

Thus, the relation between stress and strain becomes: 



A*. o> 



1 
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U-21 ) 



t'= ^ * 2 ^) + (-a,*-i.v,lj£j 

^ + n co ^ J 



+ VY) 1 1 



~f' r ^f Pf , oE»+ AI,f 3 •+ +■ *(im+A)I,w 



■l-Mcoy| +- ] 



or for 
(B-22) 



fo - o 

T'= & [M,£ 3 i-ZG^-l-(^r^-ZmI x )l£ 3 




Murnaghan [ij, has shown that 

(B-23) olei (E 3 +2y) = (del -TJ* = | + 21, 4-41,. +8l i 

01 fL = ( I *2t MI, *■ ?0 1 



To a second order approximation, the compression ratio is 
(B-24) 



%- (i + 2i, «• <*i,) '*• = i-r, 

neglecting terms higher than 2nd order, we get 

t = m,£s+ 26^ + ^ c »»+*) I. ^ + r\ co y 6r,^ 

(b- 25) T= ‘xr 1 E 3 +- 2 ^ t- et 3 +- z (m- x- &) L ^ 



<n 



+ n co 



+ ^ Cr *l>l 
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SUMMARY 



Third-Order Approximation of stress-strain relation (three 
dimensions ) . 

( 13 - 25 ) + t-ncoj + 46 -^ 



where “]“ 

E 3 

I, 

lx 

l, to, *1 
> 

y 

E 

Gr 



= stress tensor 
= strain tensor 
= identity tensor 

= First invariant of strain tensor 
= Second invariant of strain tensor 
= Second order elastic constants 
(elastic constants of Lame) 

= Third Order Elastic Constants 
= 

G+M-**') i -ay 

= Poisson's ratio — — _ . 

= Young's Modulus = G-(3?vf 26r) = 2& (f-+y) 

■>vt-Gr 

= g = shear modulus 

ao+y) 
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APPENDIX C 



DERIVATION OF CLASSICAL EQUILIBRIUM EQUATIONS 
I. Balance of Forces 

Considering the elemental plate shown in figures C-1 and 
C-2, we have for the x and y component forces 
x component 

(c-i ) (N* +• oU)clu “ Nx clij +■ ( ^ cAij.) c/x 

- dxd<j. 



= dx 

Sx J 



= +. ^_n^a. +X = o 

d * & Lj 



where X is the x body or tangential force component per unit 
area of the middle plane of the plate 

y component 

Similarly we have for the y component forces 
(C-2) *>_N^ + + Y = O 

<^X ^ Lj, 



If we assume and are small compared to dx and c/ u » 

£ x £ 7 

we have for the projection of the N ^ forces on the Z. axis: 
(C-3) *■&*>)<* 1 - N*|* 

- ' ^ d/du. + N X ^-\o(xciu -b^H* <A\ck<± 

e* QX I dx x / 27- A J 
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or, reducing to second order terms, we have 



(C-4) N x d X cl 

0 X x 



1 



and similarly, the z component of the N forces becomes 
(C - 5) •»**, 

since = Ni^.y. from equilibrium of moments about the Z. axis 

the Z component of the Nx^ and N^x forces is 

(c ~ 6) 2 Nxh iV 

' 0X0 



1 



Thus, the membrane portion of the Z component forces become 

(C-7 ) N* kw +- +Mu ^ 



The effect of shears (3* and , as shown in figure C-3 also 

contribute to the equilibrium of Z component forces, and can be 

shown to be: ( 4 - ) dx c(q 

' 0 x ^ ij. / J 

The lateral load component is ^9 dx d<-^ and the body force com- 
ponents are a„4 

* 

Thus, the combination of Z components in the general case is: 

(c-8) P + iQx + -Y 

' &y. z <r 



C-2 



<M<2x 



BALANCE OF MOMENTS 



Considering figure C-3, by taking moments about the x axis 
and considering the right hand rule for positive moments, we 
have 




‘1 



4 . ij +■ 

+-(( 2 ^ +■ dx. - 



o 



Simplifying and retaining terms up to the second order only, 
have 



we 



(c-io) - iiil h- Qu - o 

£ X ‘y j 



Similarly, by taking moments about the y axis, we have 



(0-11) SMx , 
c>x 

Differentiating the 
2nd with respect to 



^4*- - Qx =o 

1st equation with respect to 
x and adding, we have: 



y and the 



( C— 1 2 ) 


£ a Mx 


4 - y M a x 4- 


- 


^*4 


o) (2 tA __ 




X * 






b * &<y 


^ ij. 


or since 


'V 


■ - M n 


(?■.. j 


, r r ) 




we have 


(C-1 3 ) 




_ 


4- 


— 


+• ^ 








&>X 





^Qx 



o 



But from the 2 component force equation, we have 

(C-14) = -/f -t-Nx^ a. Nu -X ^ - V" £*/ ) 

&x ^ \ ' tx 1 
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Thus we have 



( C-1 5 ) 



- - 4> fNx^v + iNrMik 

& **■ <5^* y' ^x 1 ^ 



The moments resulting from the stresses distributed on the x and 

y faces of the differential element shown in figure C-3 are 

, t/x 

( C-1 6 ) M* * j Zclz 

-Vi 



(C-1 7 ) 


m 5 = 


( t/x 

J 6uZdl 

-Vi - 1 


(C-1 8) 




r i '- 

= " J Zxnhii 

“Vi. J 


From the 


plane 


stress equations of 


we have 






(C-1 9) 


= 


g («* 


(C-20) 


e, . 


e ( <s^ - y 6\ ) 


( C-21 ) 


Y ") = 


zo+y) ■>* 
— ^ 


or, in terms of 


stresses, we have 


(C-22 ) 


6* = 


Ce^ye^) 


( C-23 ) 


6 1 = 


^(e 3+ xe,J 


(C-24) 


r n 


- O- V"xcj — 



_JL_ r*w 

2(i+y) -> 



C-4 



since 



e* = 


z — 

r* 


__ ^ &**/ 
21 S>x* 


e x = 


z - 


_ -* S 1 */ 


Vxij- " 


~ Zz 





by substituting the stress relations (C-22), (C-23) and (C-24) 
in terms of Z and the curvatures into equations (C-16), (C-17) 
and (C-18) and integrating, we have: 



(C-25) 
(C-26) 
(C-27 ) 



M, 



M, 



^ 1 d q w 






_ f . y \ 

- ^ ( d\i l z / 



Mku ^ D(i-/) -^ x ^- 

- 1 <^3 



where D is the flexural rigidity 



D = 



El 3 



\iC\-y x ) 

Differentiating equations (C-25), (C-26), and (C-27), we have 



(C-28) ^M x 

£ *"*• 






^ >> 






:0 



(C-29) ~ 



(C-30) 



Thus 
(C 



*>X ^li 



a i _ D(i-y) 






-31 ) +ZO->')£k ) 



-pv w 
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Finally, substituting equation (C— 31 ) into equation ( C— 1 5 ) we have 

7 0 



(C-32 ) V 4 W = f- 4- 

0 D 



qua xi on v ^ — ^ i / m iu ^4 uct ^ ^ ^ 

(N x + .? n* h + Mu _"V~ o)w _-v - ^ ^ 



Stress 


Function F 


defined as 




^ <$x 


= Nx 


V 




t 


yp 

o>x^ 




_ Nu 

t 


_ £lE 

t^X 




_ N x>j 

t 



¥e have 



( C-33 ) V W = ~—( ^ ^ ^ w — 2 -f 9 LE VW _ X ^ K/ _^V~ c)w 

D t> 6 
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Deformed Plate Element Shewing Forces and Slopes ( Greatly Exaggerated ) 




Figure C-3 



Sub-Element Showing Moments and Shear 
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